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ABSTRACT 

We consider some radiative corrections to the lowest order annihilation di¬ 
agram for the orthopositronium decay rate. The insertion of the renormal¬ 
ized vertex correction in the annihilation graph gives 1.6283 (a/vrj^To. We 
compute also the contribution of the square of the lowest order annihilation 
amplitude, which turns out to be 0.1702 (a/Trj^To. Finally, we obtain a term 
In a To arising from the correction to the light-light scattering block due 
to the exchange of one coulombic photon, in agreement with earlier compu¬ 
tations. 
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1 Introduction 


Measurements of orthopositronium (Ops) decay rate in the recent years 
pose a great challenge in QED due to the large discrepancy between the 
experimental result and the theoretical predictions [|l|. In fact, the two most 
accurate experimental rates ii 


= 7.0514 ± 0.0014 fis-^ and = 7.0482 ± 0.0016 /rs“^, 

deviate by 9.4 a and 6.2 a from the theoretical valueQ, whose most accurate 
estimate to order o/tt has been given in P: 


Ag^f(P^^ 37) = 7.038236 ± 0.000010 
The decay rate to leading order, 


Tn = 


a^mc^ 2 (vr^ — 9) 


= 7.21117ps 


-1 


( 1 ) 


h 9 TT 

was computed by Ore and Powell [Q and the order a corrections by several 
authors 0, |, g 0 §. 

Faced up to this difficulty, theorists have made a great effort to compute 
the next correction, (o/tt)^, of the perturbative expansion [^, [T^, 


Ag^f (P^^3 7) 


= Pn 


+ O 


1 + (-10.2866 ±0.0006) 


a 


TT 


a 


TT 


±0((a/7r)3) 


( 2 ) 


To get rid of the theoretical-experimental discrepancy, the coefficient C must 
be of order 250 ± 40. Such a large value cannot be excluded, even if it may 
appear unnatural in the framework of perturbation theory (PT). If the results 
of [0 are correct (see footnote |I|) only (P ~ 30 is required. 

The computation of C is very hard due to the large number of Feynman 
diagrams contributing to the order of PT. Some of these have been already 
calculated: the vacuum polarization type corrections to the hrst order graphs 

^Recently a new experiment Q gave the value 7.0398 ± 0.0025 ± 0.0015 (the first 
error is statistical and the second systematic), in good agreement with the theoretical 
expectation. However, an independent confirmation of this measurement is necessary 
before concluding that the Ops problem is experimental instead of theoretical. 
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were considered in 
ing block in 


rT|, the radiative corrections to the light-light scatter- 


1^], and the square of the first order amplitude in m- The 


relativistic corrections, i.e. those associated with the expansion in n/c ~ a, 
where v is the relative velocity of the e’''e“ pair in positronium, were taken 
into account up to order in ref. [|T^ . 

All the contributions to the amplitude up to second order PT may be 
written in the form: 


O' 

M = Mq + — ( Mg' + Ma + Ml ) + 
vr 

( Mab' + Mar + M 2 ) + 0 , (3) 

whereQ Mi represents the sum of all the hrst order amplitudes with the 
exceptions of the annihilation diagram (see £g. 1), denoted by and the 
subtracted binding amplitude, Mg/ |^. The second order annihilation type 
corrections are given by the subtracted binding diagram, (fig. 2-A) and 
the radiative corrections to the light-light scattering block. Mar (an example 
of which is given in fig. 2-B). M 2 denotes the remaining (non-annihilation 
type) second order amplitudes. 

In this paper we consider the second order corrections given by Ma and 
Mab' and the logarithmic enhancement produced by a coulombic one photon 
exchange in Mar- The contribution of these corrections to the decay rate 
has the form 


a 


vr 


2 Re{M*g, Mq) + |M^|^ + In a \Mq 


( 4 ) 


The previous expression must be summed over the final photon polarizations, 
averaged over the Ops spin states and integrated over the phase space of the 
three hnal photons, with the proper kinematical factors (see for example |^). 

The remaining part of the paper is organized as follows: in the second 
section we compute the contributions of the binding corrections to the lowest 
order annihilation diagram. In section three, using the known results for the 
light-light scattering tensor, we compute the contribution of the square of 
the lowest order annihilation amplitude. Finally, in the fourth section we 

^Note that here and in the rest of the paper we write explicitly the powers of a/n, 
relative to the lowest order Mg, for each amplitude. (For example, Ma is of the same 
order as Mg). On the contrary, we omit them in the text. Note also that we will not write 
any power of a/ir for the two unsubtracted binding amplitudes, Mb and Mab, since they 
contain terms of different order in a/n. 
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consider the exchange of one conlombic photon in the light-light scattering 
block of the annihilation amplitnde, hnding a logarithmic contribntion in 
agreement with earlier calcnlations B © 0- 

2 Radiative Corrections to the Annihilation Process 

It is well known p, ^ that the contribution to the amplitude for the Ops 
decay rate originated by the binding diagram contains also the lowest order 
approximation: 


OL 

Ms = Mo + -Ms'. (5) 

TT 

Therefore only the subtracted binding diagram Mg' must be included, oth¬ 
erwise Mq would be counted twice. This phenomenon occurs due to the 
presence of the conlombic part in the virtual photon propagator, which had 
been already taken into account when solving the Bethe-Salpeter equation 
for the Ops wave function. 

It is quite clear that the analogous phenomenon should come out in the 
“binding” type radiative correction to the lowest order annihilation diagram 
(see £g. 2-A), and in fact we shall show in this section that the amplitude 
Ma is contained in M^g. Therefore, we can write: 

M^s = -M^ + f-V M^s', (6) 

TT V VT/ 

where M^g' is the subtracted binding-annihilation amplitude. 

We express the amplitude as a product 


AB 


4m2 ^ 


( 7 ) 


In the previous formula the 4-vector describes the transition of the 
heavy photon to three real ones, A = (A 2 , A 3 , A 4 ) stands for the set of the 
three polarizations of the hnal photons, Aj = ±1, and is the order a 
correction to the annihilation current 4-vector of the positronium in the 
polarization state em- Explicitly: 


P ~ 47r 


d^P fd^k A^,{k-p) 
(27r)^ 7 {k — pY 


X 


Tr I 7p —P/2 + k + m 7^ 


[ {-P/2 + kY-m‘^][ {P/2 + kf - np] 


P/2 + k + m 


7< 


4 


+ Ct , ( 8 ) 
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where Ct stands for the contribution of the vertex counterterm. We use the 
notation k = and here and everywhere in this paper a term iO must 
be implicitly understood in each factor of the denominator arising from a 
propagator (only in eq. (|^ the iO will be explicitly written). In eq. (§) P 


is the Ops momentum in its rest frame: 


m ol 
2 


P = {2W, 0, 0, 0), with W ^ m — ^ , and 7 = 

The wave function of Ops relevant to our approximation is (see for 

example P, PT[]): 


= {27r) 6{po) V2m 


0 

0 


a ■ er. 


0 


<h(p), 


with <h(p) the nonrelativistic ground-state wave function 


$(^ = 00 


8 Try 


1^2 ^ ^2y 


(9) 


( 10 ) 


and the constant 0 o is the wave function at the origin, 0 o = 

The tensor obviously depends on the gauge we use. The choice of 
the gauge is subtle when dealing with bound state problems. It has been 
discussed elsewhere (see for example [^) that the Coulomb gauge is the 
most natural for calculations in positronium. However, covariant gauges are 
simpler for computing radiative corrections, and, among them, the Fried- 
Yennie (FY) gauge is the most convenient, due to its good infrared behaviour. 
We shall compute both in the FY gauge and in the Coulomb gauge. As 
expected, the result is the same in both cases, and no gauge correction term 
must be added when using the FY gauge. 


2.1 Fried-Yennie gauge. 

The FY gauge is a covariant gauge dehned by 

+ 2^. (11) 

It has good infrared properties, allowing us to work safely with a zero hcti- 
tious photon mass from the beginning. 

Following IP] we separate the trace entering the integrand of (P) into 
two pieces, one which remains non-singular at k =0 and one containing the 
contribution of the coulombic photon. To this end, we dehne 
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m 


7p 


P/2 + k + m 7 j .| 


tr^yp{k) = Tr [-P/2 + fc + 

and write: trp^p{k) = trp^p{0) + {trp^p{k) - trp^p{0)}. 

We consider first the term trpi^p{0) . The 7 matrices algebra leads to: 

0 a ■ e 


( 12 ) 


tTp^piO) = -AW^ Spod^oTr l-fp 


0 0 


+ 0{ap. (13) 


To arrive to the last expression we have used the fact that 


T^"*'(p) = 0(a'^) Qf + = 0(a‘^). (14) 

The integration in pq in (H) is trivial, using the delta function entering 
the formula (^. It remains the following integral: 


f (Pp S'lT'y 

J (^ {\p\^ + Pf 


Pk 





(l + 2 j 



{k — pY 

\k + \PY-m? 

\k- 

\PY - rP 
(1 


5) 


whose result, ^ — 3 + 0{P), can be found in [Q. Hence the contribution of 
this term to is: 


2 

am ^ 
— I - Tr 

71 


Ip 


0 

0 


a ■ er, 
0 



(16) 


Let us now consider the remaining term, trpypik) — trpyp{Q). In this 
case the integral is free from infrared singularities and we can put p = 0 
in the loop integral, introducing an error of order o?. However, ultraviolet 
divergences are present; we regulate them by using dimensional regularization 
(the analogous calculation in cut-off regularization is performed in appendix 
A). It is important to remember that some care is necessary when using 
dimensional regularization in the FY gauge. As was shown by G. Adkins 
[|T7[| , it is convenient to choose the tensor as: 

= 9,. + (17) 


where e = {A — d)/2 and d is the complex space-time dimension. 
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Using 7 matrices algebra and (|^ one can see that the only contribution 
of tr^ypik) — trp_^p{Q) to the integral (|) is given by: 


rp j r r 0 CT ■ 

Trjjpkjpkjp Q Q 

In this way we arrive to the following integral: 


(18) 


iea 

f d^k ^pk^pk-iy /2^pu{k) 

dvr J 

i7r2(27r/i)“2^ 

{k -|- \PY — 

\k-\PY-m? 


where /i is the dimensional parameter introduced in dimensional regulariza¬ 
tion. The integral ([I^) can be evaluated by standard techniques, giving 

(3D + 8-2dpo)7p, (20) 

where D = f— 7 £;-|-ln (The number = 0.57721 is the Euler constant). 
Note that current conservation implies PP = 0. Since we work in 

the positronium rest frame (P = 0), it follows = 0 and the term with 
5po in (^) can be ignored. The contribution of the remaining term to 
is then: 


za 


- (3D -h 8) Tr < 7 p 


0 a -eZr 
0 0 


( 21 ) 


Adding ([l^) and (|^) , inserting the result in (J^) and noting that the hrst 
order annihilation amplitude can be written as 



71 



0 a -e^ 
0 0 




we obtain the following contribution to M^b 


( 22 ) 



TT 


1 + - 
TT 



3D + 8 Y 
4 j. 


(23) 


The contribution of the vertex counterterm will cancel the divergence ap¬ 
pearing in (pSl). Due to the Ward identity, this counterterm can be obtained 
from the self-energy correction to the electron propagator. 

To order a the mass operator has the form: 


7 



















(24) 


m 


_c^ r ^ 7^ (/ - fc + m) 7 ^ 

Att J ztt^ /c^ [(Z — ky — m^] 


On mass-shell renormalization conditions imply 


iZ 


p — mo — Tj{p) p — m 


p 


m 


(25) 


where mo is the bare electron mass. From (|2^) it is possible to obtain the 
electron wave fnnction renormalization constant, which in the FY gange is 

inn: 


Zfy = I - ^ (3D + 4) . (26) 

47r 

It follows that the contribntion to the amplitnde of the vertex connterterm 
is: 


/ay 3D + 4 
vvr/ 4 


M^. 


(27) 


Finally, snmming up all the contributions considered here (the terms 
trij,up{0) and tr^^pik) — tr^yp(O), and the counterterm insertion) we get: 


= 



(28) 


2.2 Coulomb gauge. 

We will show now that the same result is obtained working in the Coulomb 
gauge. This gauge is obtained by making the substitution: 

(fc _p )2 Gp^k-p), (29) 

where Gpu{,q) is the Coulomb propagator: 


Goq{(i) — — , Goi{q) — Gio{q) — 0 , Gij{q) — — iSij — ^ 


q 


q 


q 


(30) 

Again, using 7 matrices algebra and ([T^), we rewrite the trace in the 
numerator of (||) as 


tr^upik) = -PpPyTr [A!^'"\p)p] + Tr f T'’"^(p) 7 ^ fc 7 ^ ^ 7 ^ 















+ Tr 




Py'lpk'l^,) 


(31) 


It is easy to see that the last term vanishes after contracting the Lorentz 
indices with the ones of the photon propagator. The contribntion of the hrst 
term to is 


— (IGtt) Tr 


7p 


0 a -eZr 
0 0 


where Jq is defined by: 


(32) 


/n = 


Pp 8 Try 


(2vr)3 + ^2)2 

d'^k 


X 


(dvr) 


{k 


p) [(-PI2 + kf - np] |(P/2 + kf - np] 


(33) 


This integral has been studied in |21|. Its result is 


/o 


i 1 / TT 

(dvr) 2 m? V o. 


(34) 


Using this result and considering (J^) and (p^), we find that the contribution 
of the hrst term of (^T]) to the amplitude is 


(1 - 2a/7r) (35) 

which is the total result obtained in the FY gauge. 

Now we will show that the contribution of the second term in (|3TD exactly 
cancels against the contribution of the Coulomb gauge vertex counterterm. 
Hence, the result in the Coulomb gauge is the same as in the FY gauge. 

The second term of (pTD gives raise to an UV divergence. Again, we 
choose dimensional regularization to give a meaning to the loop integral. We 
work in d = 2a; dimensions, with one temporal and 2a; — 1 spatial dimensions. 

As in the case of the FY gauge, there is no infrared problem for this term, 
and we can put p = 0 in the photon propagator of the integral (^. We have 


d 2 ^fc 


'Juk'jpk'jp 


(27r/i)2^ [{-P/2 + kf - m2] [(P/2 + kf 




(36) 


The formulae for integrals of non-covariant functions in this dimensional 

After standard computations 


regularization prescription can be found in 
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and taking into acconnt the fact that the terms with the Lorentz index p = 0 
do not contribnte to the amplitnde, we obtain: 


i 

(dvr) 2 



for the case of the temporal propagator, Goo, and 


(37) 


i 

(dvr) 2 



(38) 


for the contribntion of the spatial components Gij of the photon propagator. 
The total resnlt is therefore 


£-D(-iej^). (39) 

The contribntion of the vertex connterterm is given by 6Zi (—i ejp) . 
Since the vertex connterterm to order a in the Coulomb gauge is 

SZ, = -£d, (dO) 

we see that the connterterm exactly cancels the contribution of the second 
term in (^), as we wanted to show, and then, the amplitude is given by (p^) 
also in the Coulomb gauge. 

2.3 Conclusion. 

As claimed at the beginning of this section, we must consider only the 
subtracted amplitude, 

o. /ci \ ^ 

-= -2 - . (dl) 

TT VTT/ 

Hence the Ops decay rate receives a contribution given by the integral of 

i S 4 E 2 fle (m^« • M<”'') (42) 

over the phase space of the three hnal photons. Since Mab' is proportional 
to the lowest order annihilation amplitude, the contribution to the width 
is proportional to the lowest order annihilation width, whose value can be 
found in p]: T^ = —0.81d05 (ck/tt) Tq . 

Therefore, the renormalized vertex correction to the annihilation amplitude 
turns out to be: 


10 




( 43 ) 


r^B. = -2-rB = 1.6281 ^r„ 

TT Tl^ 


Note that, differently from the lowest order annihilation case, it contributes 
positively in the direction of reducing the theoretical-experimental discrep¬ 
ancy. Its numerical value, however, is too small to make a signihcant progress. 


3 Lowest Order Annihilation Diagram. 

Now we consider the lowest order annihilation matrix element (fig. 1): 


a/J 


TTCK 


d^p 




(44) 


4m^ J (27r)‘^ 

The integral over p is trivially performed. The next stage is to square the re¬ 
sulting matrix element and make the average over the Ops polarization states 
[§]. After simple calculations, similar to those performed in the paper of one 
of us 0. — where the production of tree gluonic jets in electron-positron 
colliding beams was considered—, we obtain the following contribution to 
the width^: 


T .2 = 


->"^15^ J rf=’i^6(E^.-2) E,G“G™^ (45) 

To arrive to (^51) , we have expressed the phase space volume of the three 
hnal photons as: 


d^ko d^k^ d?k. 


U)2 Co's Ct'4 


6‘^{P—k2—k3—k4) = STT^m^ y (ii/2 (iz/3 dz/4(5 (2—1/2—z/3—1/4), 

(46) 


where Ui = ^ . 

'' m m 

At this point we use the results of papers [jT^, namely: 


= 2®a^ [i?(234)+ i?(324)+ i?(423)] ; (47) 


^Remember that the vector is space-like, as explained before. Therefore, is 
positive, in spite of what at first sight might seem due to the minus sign in the r.h.s. of 

il)- 
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2 


/2(234) = i?(243) = 
i £i"i+(234) ' + 


£+++(234) 


+ 


^2 


Z/3l/4a2 


£-++(324) 


+ 


1 


^S+(234) + £. 


S+(243) 


+ 


( 2 - 3(24 
1/2 02 


-£S+(234) - lfS+(243)'' 

O 3 04 


( 48 ) 


where a* = 1 — z/j. The rather cumbersome functions T, whose arguments 
are the z/j, were calculated in the paper of Costantini et al. [^. Their 
expressions are explicitly written in appendix B. After numerical integration 
over the phase space, we hnd: 


T 



(49) 


/ 6{2 - EiUi) [ R{234) + i?(324) + i?(423) ] 

32(7r2 - 9) 


0.17021(10). (50) 


We would like to stress that the contribution from muon (and hadrons) 
as fermions in the loop is negligible (of order {me/m^Y) [^. 

In conclusion, the total correction hnd here, adding (^) to (H), is 
1.7983 (a/7r)^ro. Manifestly, it is still far from solving the discrepancy be¬ 
tween the modern theoretical and experimental results. If the Ops problem 
is to be solved by this kind of perturbation theory, larger contributions to 
the width must be searched another class of diagrams. 


4 Coulombic Exchange in the Light-Light Scattering Block. 

In this section we shall obtain a known logarithmically enhanced contri¬ 
bution arising from radiative corrections to the light-light scattering block. 
The whole set of these radiative corrections has been already computed by 
Adkins and Lymberopoulos ||^. We present here a simple computation of 
the term of order In a, which arises from the diagram displayed in £g. 
2-B. 

The contribution to the amplitude of the radiative corrections to the 
light-light scattering block can be represented by: 
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M'» = -le J {>P<”’(p) 7 ,} ^G<»f{k,,k,,h) , (51) 

where the subscript (R) in means the order a radiative corrections. 

We are interested in the contributions to of the loops of the form 
represented by £g. 2-6.^ It is given by 

ie^ r dH r d^k (/c — Z) 

(dvr)^ J (27r)^ J {I — ky 


Tr 


■lu 


— \P + l — m —\P + k — m 


qw 

OQ/1 


234 I 


7/. y- 


P + k — m ip + I — m 


■7p 


(52) 

In the previous formula — 1) is chosen in the FY gauge (pT]) and 

describes the annihilation of a pair e"*" e~ at rest to three photons, namely: 




Te\, 


jP + k + k4 — m 


1 ■ ^As 


ip + k — ko — m 


7 ■ cas . (53) 


It is worthwhile to underline that there are two regions in the loop mo¬ 
menta space from which the integral (^) receives the main contributions: 
one, where the fermion momenta are far off mass-shell, which correspond 
to Z, fc ~ m, and another one with fermion momenta almost on mass-shell, 
l,k ^ am. It is this last region which originates the logarithmic term. In the 
analysis of this region we can neglect k and Z in the numerators of fermion 
propagators (after rationalizing them) and in Then, using 7 matrices 

algebra, we can rewrite the trace in (^) in the following way : 


- P,P, Tr { Oill (1 + 70) 7p (1 - 7o) } + O(a'), (54) 

where we have used the fact that P^ — 4m^ = O(a^). The integration in k is 
now identical to that of ([I^). We use the result [§]: 


d'^k 




-m?l\oQ{k - 1) 


{k — iy {—^P + ky—m^ Gp + ky — 


TTm f 

^ arctan-3. 


7 


(55) 


^Strictly speaking, Fig. 2-B represents only a typical loop. There are five more graphs 
of the same kind corresponding to permutations of the three final photon states. 
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The integral in I can be performed in two steps. First, we integrate /q 
using the method of residues: 


r d/o 

1 


i 1 

rs_; 

/ 

(—+ /)2 — + iO 

'(|P + /)2 - m2 + i0' 

4m fs + ^2 


(56) 

The remaining integral over the spatial 3-momentum I is logarithmically 
divergent. However, we do not worry here for the ultraviolet sector, which 
is outside the integration region we are considering and gives no logarithmic 
contribution in a. The considerations made at the beginning of this section 
permit us to replace the upper limit by m in order to extract the logarithmic 
term which arises from the infrared behaviour: 


f d^l i 1 rmr 

J (27r)3 4^|/|2 + y ^ 


arctan ■ 


7 



(57) 


where we have approximated arctan Y to |. 

Collecting all the factors, we have as a contribution to : 

7^ E Tr{ (1 + 7o) 7p (1 - 7o) } , (58) 

' ' ctGSs 

where cr is a permutation of (234). The amplitude receives a contribution: 


AA/r/> = a"loga7y'7^rr{4."">(p)77 X 

y rr{0«„,(l+7„)7,(l-7o)} , (59) 

ctGSs 

and the corresponding correction to the width is the integral over the three 
hnal photon phase space of 

A E 42 * K ‘"') ■ (60) 

0^0.^ 

The lowest order amplitude, Mq, can be written in our notation as 

AC'’ = E / T^rr {4-'“>(p)OtU} ■ (61) 

0-GS3 '' 
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Using techniques similar to those described in P , it is possible to rewrite 
the product of the traces in such a way that the leading width Fq appears 
explicitly. We hnd 


AF.r = In — Fn . 


a 


(62) 


This result is in agreement with previous calculations 
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A Appendix A. 

In this appendix we discuss the contribution to of the term tr^upik) — 
and of the vertex counterterm, using as a regulator a sharp cutoff A 
in momentum space. Of course the physical prediction coincides with that 
obtained in section 2 in dimensional regularization. 

As discussed in section 2, the part of the integral (|]) given by trpj,p{k) — 
trp,up{Q) has no infrared singularity and we can put the internal positronium 
momentum p equal to zero in the loop integral, making only an error of order 
a^. Hence the integral in p can be trivially performed; it remains the loop 
integral. Introducing Feynman parameters, making simple integrations inQ 


— le 


a 

dvr 




'jpk'jpk'j^ A 

flU {k) 


(k + lpy 




{k-\py 




(A.l) 


and ignoring a term proportional to Pp, which does not contribute to the am¬ 
plitude Mab due to the gauge invariance of the light-light scattering tensor, 
we get 


2 


Ati 


7p 


Therefore the contribution of this term to the amplitude is : 


(A.2) 


a 

dvr 


3 In 


A2 


+ 






(A.3) 


The computation of the self-energy —to get the vertex counterterm— is 
more subtle. The mass operator to lowest order is: 


m 


.a r d^k'y^(l-k + m)-ipAp^{k) 

^ dvr J iTT^ {k^ - A2) [(/ - ky - m?] ’ 


(A.d) 


where A is the hctitious photon mass. Note that in the FY gauge the numer¬ 
ator does not contain any term which produces linear divergences, therefore 
the Feynman trick to join the denominators and the subsequent shift of the 
loop momentum can be used, giving the result: 

®In this case Apv{k) is given by (0)- 
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E(p) = 


— la 


dx 


/ A2 \ 

(3m — p{l + x)) In — — 1 — 


271 Jo 

^ [{1 — x)a^ + Ax^{l — 


2a? 


(A.5) 


with 


A2 


p 


m = m^ (x^ + z/(l — x) -\r x(l — x)p ), u = —- , P = 1-5- • (A.6) 






The electron wave function renormalization constant, in the on-mass-shell 
renormalization scheme, is defined as 


Z = 1 + 


d^P) 

dp 


(A.7) 


Simple algebra leads to the following expression for Zpy'- 


a 


1+ — / dx 
277 Jo 


where 


A2 


9 


(1 + x) In —- — 2 Ina: — 1 — -(1 — x) -\- 2J{x, u, p) 




>0 

(A.8) 


J{x,u,p) = 


a; (1 — x)^ (2 — x) (z/ + px) 


[x^ + z/(l — x) + x(l — x)p]^ 

The quantity Jq Jdx depends on the way in which z/ and p tend to zero: 


(A.9) 


Jdx = 

1 

2 

p « V - 

-0, 


Jdx = 

2 

V << p - 

^ 0 . 

(A. 10) 


Y. Tomozawa il showed that the right result appears in the limit v « 
p —0. The resulting expression is then: 


a 


A' 


= 1 - — I 3 In — - - 


dvr 


m^ 


(A.ll) 
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This implies that the contribution of the vertex counterterm insertion is: 


a 

An 


2 


M, 


(A.12) 


To have the total expression of Mab we only need now the contribution 
of the term tr^yp{Q). It does not contain any divergence, hence it can be 
recovered directly by eq. (0) . We have: 


a 


a 


^ ^ n) U W ■ 

Collecting the results of (|A.13|) , (|A.3|) and ( |A.12|) , we get the same as in 
dimensional regularization (p^). 


B Appendix B. 


We write here the explicit expressions for the quantity i?(234) entering the 
tensor G. Accommodation from the results of paper to the annihilation 
channel was done in the paper ||19|| . 

For the case e = m (notation of El) we obtain: 


4++(234) = 

20304 


1^3 


+ ^(^ + ^-as) (B(as)-B(l)) + 
V O2 U 3 / 


1^3 V O2 
2 1 


203O4 (— + —'j (-8(04) — -8(1)) + 

V 02 z/4 / 


.02 z/4. 

'^Os - 04 - (7^(03) + T(o4) - T(l) - 1(03,04)) 


02 


02 


^T(o2) + - (T(o2) - T(o4)) - T(o3) + 

(2-4 (22 


204 


03 


,^3 


Vn 


1 

02 


Z /2 (02 — O4) 


O2O4 


/(O2, O3) 


O4 
Z^203 


(T(o 3) - T(l)) + 


O3 


7(02,04 ) - -(T(o 4) - T(l)) + 

02 04 z/4 


2 _ ^ ^ ) 7(03,04) ; 

04 02 


(B.l) 


fti+(234) = 
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03 [— — —^ (^(^ 2 ) + T(as) + T(a4) — T(l)) — 

\ 04 02 / 

— 1(02,03) + —1(03,04); (B.2) 

04 02 


£S+(234) = 


4o3 


2oq 


4o4 


204 


^ ( 5 (^ 3 ) _ 5 ( 1 )) + ^ ^ (5(04) - 5 ( 1 )) 

02/23/ ' - ' 

/ 40304 ^ 

V «2 02 

1 1 


O2 /24 

(^(03) + ^(04) — ^( 1 ) ~ -^(03504)) — 


- + - + - ) 7^(02) 

(I 2 Cts <^4 / 

(^ + 4) T(04) + ( 

V/2403 02/ V 


Vz 


T(o3) 


Z/4O3 
O4 1 
O3O2 
Z/2 + 1 
(24CI3 


+ 


Z/3 

Cl2^4 


1(02,03) + 


02 ^ 

Z/304 02, 

--- + -)T(1) + 

O 3 O 4 z /3 Z /4 02/ 

O3 + 1 ^ 


Z/2 


Z/4 


O4O2 


O2O3 


7(02,04) + 


+ —1 7(03, 04) ; 
02/ 


d^i+(234) = 

—2 — — + — + —^ (7"(o2) + 7"(o3) + T(o 4) — T(l)) + 

V 02 03 04 / 

f— + -1 7(02,03) + f^ 7(02,04) + 

\04 O2O3/ \03 O2O4/ 

f— + -1 7(03,04) . 

\02 O4O3/ 

In the previous formulas we used the notations: 


di 

B(z) 

T(z) 
I {as, 02 ) 

F(o,7) 


1 - z/j, z/2 + z/3 + z/4 = 2 , 


— 1 + Y-1 arcsin , 

— (arcsin V^)^, 

F(o3, 7) + 7^(02,7) - 7^(1,7), 


'•1 dx 


/o 7^ — 


In 


1 — 0(1 — x^) 


(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 
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Figure 1: The annihilation graph. 



(A) (B) 


Figure 2: Two different kinds of corrections to the annihilation graph. (A) 
is the vertex correction, (B) represents the insertion of a photon into the 
light-light scattering block, which generates a logarithmically enhanced con¬ 
tribution. The direction of the fermions in the final loops is clockwise for all 
the two diagrams. 
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